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Superconductivity Driven by the Interband Coulomb Interaction and 
Implications for the Superconducting Mechanism of MgB2 
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Superconducting mechanism mediated by interband exchange Coulomb repulsion is examined 
in an extended two-band Hubbard models with a wide band crossing the Fermi level and co- 
existing with a narrower band located at moderately lower energy. We apply newly developed 
path-integral renormalization group method to reliably calculate pairing correlations. The cor- 
relation shows marked enhancement at moderate amplitudes of the exchange Coulomb repulsion 
taken smaller than the on-site repulsion for the narrower band. The pairing symmetry is s-wave 
while it has unconventional phases with the opposite sign between the order parameters on the 
two bands, in agreement with the mean-field prediction. Since the band structure of recently 
discovered superconductor MgB2 shares basic similarities with our model, we propose that the 
present results provide a relevant clue for the understanding of the superconducting mechanism 
in MgB2 as well as in this class of multi-band materials with good metallic conduction in the 
normal state. 
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Mechanisms of superconductivity have been a long- 
sort issue since the discovery in 1911 but are still un- 
der active studies with aiming at novel mechanisms and 
searching for materials with higher transition tempera- 
tures. In particular, mechanisms driven by the electron- 
electron interaction in general are widely and intensively 
studied issue since the discovery of the high-Tc cuprate 
superconductors JlP although there exist not many su- 
perconducting compounds that have been established as 
driven by the electron-electron Coulomb interaction. The- 
oretical studies have intensively focussed on the single- 
band correlated electron models. In this letter we discuss 
a different possibility of a high-temperature supercon- 
ductor by focussing on a wide metallic band coupled by 
a moderate amplitude of exchange Coulomb interaction 
to the narrower coexisting band or localized orbitals. 

When two different bands or orbitals are located near 
the Fermi level in metals, several completely different 
physics may appear. One is the Kondo effect and the 
heavy fermion behavior where the hybridization and an- 
tiferromagnetic spin exchange interaction between local- 
ized and itinerant bands play a role for the spin sin- 
glet formation and heavily renormalized effective mass 
of the quasiparticles in lattice. On the other hand, if the 
hybridization is vanishing or small between the two or- 
bitals, while the exchange Coulomb interaction becomes 
relatively important, the Hund's rule works efficiently 
and the atomic high-spin state generically appears, if 
the level difference between the two atomic orbitals, F is 
relatively small. With this background, the double ex- 
change mechanism works in realizing the ferromagnetic 



ground state in metals. The high-spin state appears es- 
pecially when the ratio between the onsite Coulomb in- 
teraction for the lower-level orbital, U, and F, namely 
U/r is large. If the interorbital on-site repulsion V is 
taken into account, this rough estimate of the ratio is 
replaced with {U - V)/r. 

However, if {U—V)/T becomes relatively small or com- 
parable, the low-spin state may be stabilized and the 
Hund's rule coupling becomes irrelevant. Near this sta- 
bilization point, {U ~ V)/T < 1, another term of the 
exchange Coulomb interaction between the lower and 
higher orbitals may work for the stabilized low-spin ground 
state. This term of the exchange transfers the up and 
down spin electrons on the same orbital(band) to the 
other one in pair. In contrast to the ferromagnetism 
driven by the Hund's rule term, we show that this pair 
transfer part of the exchange interaction strongly fa- 
vors the singlet superconducting order in an appropriate 
range of the parameter space. |-, 

In 1963, using the BCS approximation, Kondco con- 
sidered the pairing mechanism assisted by the interband 
exchange Coulomb interaction between s and d bands 
for the mechanisDx.of the superconductivity of transition 
metals. YamajiclB also considered the similar mecha- 
nism in a specific condition of two-band Hubbard mod- 
els. A similar interband mechanism was alsQ.examined 
in strongly correlated models by the authors.© 

Here, we reexamine the mechanism driven by the ex- 
change Coulomb interaction and show by a newly devel- 
oped numerical technique that this mechanism works for 
simple metal compounds mixtured with elements with 
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a localized orbital (or relatively narrow band) located 
moderately below the Fermi level. Its implication for the 
mechanism of the superconducting tcansition in MgB2 
recently discovered by Akimitsu et al.EP above tempera- 
ture 39K is an important subject of discussions in this 
letter. We discuss that this mechanism may drive or 
substantially enhance the superconducting transition of 
MgB2. 

Our Hamiltonian 

n^Hi+H2+H3 + H4 (1) 

consists of the conduction band Tii with the wide and 
narrow bands or localized level part with the on-site 
Coulomb repulsion terms Ti.2 and Tis and with their in- 
terorbital Coulomb interaction 7^4 = Ti^E-'r'H.iP+'Hiv- 

Hi = - ^ i^^J(cLc^i^ +H.c) 

-fiNi -{(1 + T)N2 (2) 

7^2 = Ui^nu^niii (3) 

i 

^3 = C^2^?^2JT'^2Ji (4) 

j 

UiE = ^J^^ c\-^Cii„'cli„iC2icr (5) 

'Hap = J ^(1 - 6„„i){c\^^,c\^^C2iaC2ia' + H.c.) (6) 
Hiv = ^ X! "li^sj (7) 

3 

N = Ni+N2, (8) 
N,=Y,nua, (9) 

N2=Y,^2la, (10) 
la 

In Eq.(l-lO), the annihilation (creation) operator of the 
orbital l{l = 1, 2) at the site i with the spin a are repre- 
sented by ciia{cj^cr)- In this letter, we assume that Ui is 
small by considering a wide metallic band 1 while t2ij is 
smaller than tnj . 

For the region characterized roughly by t/2 — F < F— J, 
the atomic ground state of the orbital 2 is fully filled sin- 
glet. However, when J becomes not too smaller than 
r, the singlet pair can transfer to the conduction band 
1 through the term (6) and it generates a charge fluc- 



tuation with a singlet cloud around a site of the local 
orbital 2. In this letter we show that the overlap of such 
singlet clouds formed around the neighboring atomic or- 
bitals 2 may drive the delocalization of the singlet pair 
and lead to coherent ground state with the supercon- 
ducting symmetry breaking. Such formation of itinerant 
singlet pairs survives, as we show later, when the orbitals 
2 form a band through their mutual transfer. If we take 
the region with the low-spin ground state, the terms (5) 
and (7) do not play relevant roles and we neglect these 
terms in the following discussion for simplicity. 

Here we first neglect the dispersion of the lower-level 
band, t2 in (2) and the Coulomb repulsion Ui within the 
conduction band in (3) for simplicity and discuss their 
effects later. When we take the mean field approximation 
with the uniform superconducting order parameter 

-l = ^E(^HT^Hi) (11) 

i 

-2-^E(4T4i), (12) 

i 

the selfconsistent equation for the order parameters ob- 
tained from the decoupling of the term (6) for t2ij = 
is given by 

2Si = S2 J/^S^+n (13) 
2S2 = -U2E2/^El+T^ 

where the density of states of the conduction band 1 
is Di{e) with the dispersion e{k) = j tnj exp(ifc(ri — 
rj)) — j-i and Wi is the bandwidth of the band 1 in propor- 
tion to 1/ti. When the density of states is replaced with 
the rectangular constant one for the conduction band for 
simplicity, namely, Di{e) — 1/Wi, the selfconsistent so- 
lution at temperature T = in the weak coupling BCS 
limit is given by 

5i ~ W^iexp[-^(C/2 + F/a:)]. (15) 

Here x denotes the number of the orbital 2 in a unit cell 
and if the orbital 2 forms a band through nonzero t2ij, 
basically x may be replaced with the integrated density 
of states of this narrower band with the width W2 ■ 

An important point of this pairing mechanism is that 
the prefactor in Eq.(15) is determined from the bare 
bandwidth of the wide band 1. In addition, the argu- 
ment of the exponential function may have order unity 
when we may assume that J, U, F, and Wi would have 
comparable values. This may necessarily requires the 
consideration in a strong coupling region. However, here 
we took a BCS weak-coupling approach as a starting 
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point and by assuming that Wi is the largest parameter 
among the energy scales. The strong coupling correction 
is left for further studies. 

Even more important point of this mechanism is that 
the pairing order parameters on the bands 1 and 2 have 
the opposite sign, namely J2i,ji'^iia(^ii-a(^2j-aC2ja) < 
because the pairing is mediated by repulsive interband 
exchange interaction. This anisotropy may enable us to 
distinguish the present pairing mechanism from that of 
the conventional phonon-mediated BCS theory. 

This weak coupling mean-field picture is complemented 
by numerical study of the Hamiltonian (1) in the follow- 
ing. We employ recentlyJeveloped path-integral renor- 
malization group methocfl^ to perform a controlled calcu- 
lation in the ground state. Since our basic picture applies 
to any type of lattice structure, we first employ a numer- 
ically feasible square lattice for simplicity and introduce 
orbital 2 by replacing a fraction of the lattice points (de- 
noted by B in Fig. 1) regularly and assign orbtal 1 to 
the rest of the lattice sites (denoted by A in Fig. 1). We 
have introduced the B sites periodically in x and y di- 
rections with the period Ib- The wide band 1 is formed 
from the nearest neighbor transfers ti between the sites 
A while the narrower band is formed from the nearest 
neighbor transfer t2 between the sites B. The unit cell of 
this lattice becomes Ib >^h and it causes a folding of the 
Brillouin zone with multi-band structure instead of the 
band 1. However it does not alter our basic argument 
especially near the Fermi level. The transfer between 
A and B sites is switched off. We also introduce an on- 
site Coulomb repulsion U2 for the sites B (on the band 2) 
while the onsitc Coulomb interaction for the free-electron 
like band 1, Ui is taken small. The exchange Coulomb 
interaction J given by Eq.(6) is introduced between near- 
est neighbor A and B sites. The atomic level difference 
between the sites A and B is F. 

We have calculated the pairing correlation between the 
A and B sites as well as other combinations of the pairing 
between two A sites and between two B sites. Here, the 
pairing correlation is defined by 

We fix T/h = C/2A1 = 3,t2/ti = 0.3 and Ui/h = 0.2 
where the largest effect of J is expected in the region 
of comparable U2 and F as we discussed above. In a 
realistic situation, J has to be smaller than U2 and we 
show J dependence of the pairing correlation in Fig. 2. 

Figure 2 shows the amplitude of the pairing correla- 
tion between a site on the A sublattice (with the co- 
ordinate (0,0)) and all the other sites including 
A and B for a 6 by 6 square lattice with the periodic 
boundary condition, Ib = 3 and the electron density 
n = 0.5 per site. The distance between (0,0) and (i,j) 
is given by r. The interband pairing correlations P12 
has mostly negative values as we expect. The over- 
all correlation clearly shows marked enhancement for J 



roughly larger than C/2/2. The enhancement is seen in a 
comparison with a practically nonintcracting case with 
Ui/ti = U2/ti = Ji/ti = 0.01 with the same band pa- 
rameters 12/ ti — 0.3 and T/ti = 3 as before. This im- 
plies that the pairing may be driven in a realistic range 
of the parameters when a wide band is coupled by the 
exchange Coulomb interaction to a narrower band and 
the hybridization between these two bands is vanishing. 

The result also shows that the short-ranged part of 
the pairing correlation becomes larger simply with the 
increase in J while the longer-ranged part has a maxi- 
mum at an intermediate value of J below U2- This is a 
natural consequence because the strongly coupled singlet 
is tightly formed at larger J. It bears more local char- 
acter and an effective hopping amplitude of the singlet 
becomes smaller when J becomes too large as is seen for 
an unrealistic choice of J/ti = 5 > U2/ J ■ 

Although the size scaling is difficult for the moment, 
similar enhancements are also seen in larger lattice sizes. 
When the electron correlations become small, the scale of 
the enhancements becomes small accordingly while the 
enhancement itself seems to universally exist if the ratios 
of J7i, J and F are retained. The overall feature of the 
enhancement is in qualitative agreement with that ex- 
pected from the parameter dependence of the mean-field 
result Eq.(15). 

We next discuss the implication of this result for the 
mechanism of superconductivity discovered in MgB2. The 
crystal of MgB2 has so called AlB2-type hexagonal struc- 
ture with an alternate stacking-of Mg and B honeycomb 
layers. The band structureQ'ElE^P suggests that degener- 
ate bands with two-dimensional anisotropy are formed 
from B 2pcr band and its bonding part crosses the Fermi 
level near the top of the band. In the band structure 
calculations, the center of the B 2pc bonding band is 1-2 
eV below the Fermi level and two cylindrical hole Fermi 
surface appears around the axis through F to A points 
due to the bandwidth of the order of 4 eV. A rather 
separate free-electron like band also crosses the Fermi 
level, which is formed from the hybridization of Mg sp 
and B pz bands and has a bandwidth larger than 10 
eV. The Mg sp band is expected to have a rather small 
electron correlation while B p„ band may have a moder- 
ate correlation due to dominant p character confined in 
the two-dimensional honeycomb layer. Because of their 
symmetries, these two bands do not seem to be apprecia- 
bly hybridized and are clearly separated near the Fermi 
level. This band structure is essentially similar to our 
model with the band 1 (corresponding to Mgsp-Bp^) and 
the band 2 (Q-pa). Although the amplitude of the elec- 
tron correlation is not well known in this compound, the 
choice of the parameters in our calculations in compar- 
ison with plausible values of parameters for MgB2 does 
not seem to be unrealistic and provides a good start- 
ing point in the present stage. Our calculated result 
shows that the pairing enhancement becomes optimized 
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for comparable U and F while it is suppressed, for ex- 
ample, for too large F. This is consistent with the fact 
that AIB2 does not show superconductivity because the 
B Per level is deep in AIB2. On the other hand, if F be- 
comes too small, the Hund's rule coupling starts working 
and the ferromagnetic state becomes more stable. In a 
series of the binary compounds of boronides, we expect 
that the ferromagnetic and superconducting phase may 
appear rather close each other with competitions. 

In this compound, phonon modes of Boron oscillation 
may have relatively high frequencies and may take a 
strong electron-phonon coupling which prn^f support the 
phonon-mediated conventional pairing.EilEJ Our pro- 
posal gives an alternative possibility or may play a com- 
plementary role. In comparison with the phonon mecha- 
nism, our mechanism contains a specific prediction where 
the phases of the order parameter on the B band and 
B Pz bands hybridized with Mg sp should have the op- 
posite sign. This is a crucial test for the role of the 
interband Coulomb exchange as a driving force of the 
superconductivity. This may be tested in a carefully de- 
signed junction to pick up the opposite phase. Although 
the pairing, symmetry is s-wave as supported also exper- 
imentallyta' , our theory suggests that it is not a simple 
BCS isotropic pairing. 

In summary, we have analyzed a superconducting mech- 
anism mediated by the interband Coulomb exchange re- 
pulsion. The numerical simulation results suggest en- 
hancements of pairing for the combination of wide and 
narrower bands crossing the Fermi level without mutual 
hybridization but with a moderate amplitude of the ex- 
change Coulomb interaction. We have proposed that this 
mechanism may be relevant to the occurrence of super- 
conductivity recently discovered in MgB2. 
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Fig. 1. The lattice structure employed in our numerical simula- 
tion. 
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Fig. 2. The absolute amplitude of the pairing correlation \P\ 
as a function of distance for an extended two-band Hubbard 
model. The choices of J are J = 1,2,2.5 and 5. These re- 
sults are compared with the practically noninteracting system 
with Ui/ti = C/2/ti = J/ti = 0.01 The pairing correlation 
shows marked enhancement in a range of the interband exchange 
Coulomb interaction J. See the text for details. 
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